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Although it is widely accepted that most galaxies have supermassive black 
holes (SMBHs) at their centers 1 ” 3 , concrete proof has proved elusive. Sagit¬ 
tarius A* (Sgr A*) 4 , an extremely compact radio source at the center of our 
Galaxy, is the best candidate for proof 5 ”', because it is the closest. Previous 
Very Long Baseline Interferometry (VLBI) observations (at 7mm) have de¬ 
tected that Sgr A* is 2 astronomical unit (AU) in size 8 , but this is still larger 
than the ” shadow” (a remarkably dim inner region encircled by a bright ring) 
arising from general relativistic effects near the event horizon 9 . Moreover, 
the measured size is wavelength dependent 10 . Here we report a radio image 
of Sgr A* at a wavelength of 3.5mm, demonstrating that its size is ~1 AU. 
When combined with the lower limit on its mass 11 , the lower limit on the 
mass density is 6.5 xlO 21 Msun pc” 3 , which provides the most stringent evi¬ 
dence to date that Sgr A* is an SMBH. The power-law relationship between 
wavelength and intrinsic size (size oc wavelength 109 ), explicitly rules out ex¬ 
planations other than those emission models with stratified structure, which 
predict a smaller emitting region observed at a shorter radio wavelength. 
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Supplementary Figure 1 — Measured (FWHM) angular size of Sgr A* vs. observing 

wavelength. Blue and green lines represent the best-fit scattering relations of 1.39A 2 and 
0.69A 2 along the major and minor axes, respectively. Black data points (open and filled 
circles) are size measurements from the quasi-simultaneous observations in February 1997 
with the VLBA plus one antenna of the VLA. Red data points are from the best VLB1 
observations ever made, in terms of the recording rate (thus the sensitivity) and the weather 
condition, at 7 mm (open and filled diamond) and 3.5 mm (open and filled square), respec¬ 
tively. Here, open symbols are for major axis sizes and filled symbols are for minor axis sizes. 
la error bars are plotted. See Supplementary Information and Table 1 for details. 
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Details of Model Fitting Procedure for the Source Structure Determination 

An accurate determination of the apparent size of Sgr A* from the conventional VLBI 
imaging technique has suffered from calibration uncertainties. The southerly declination 
of Sgr A* (—30°) and northern latitudes for most of the existing VLBI antennas result in 
poor north-south resolution, which happens to be along the minor axis direction of the 
scattering structure. Furthermore, atmospheric effects, especially at shorter wavelengths, 
severely affect the amplitude calibration which is crucial to the determination of the 
structure of Sgr A*. 

To minimize the calibration errors and thus to improve the accuracy of the apparent 
source structure measurements, we have developed a model fitting method by implicitly 
using the amplitude closure relation . The major difference between this method and the 
widely used self-calibration technique for VLBI imaging lies in the way of obtaining the 
model of source brightness distribution. Both algorithms converge to consistent results 
with high signal-to-noise (SNR) data. The biggest drawback of the traditional imaging 
process is the non-uniqueness of the final image. This becomes more severe when weaker 
detections are involved, as in the case for the VLBI observations of Sgr A* where the A 2 - 
dependent scattering image is always resolved out on the long to intermediate baselines, 
depending on the observing wavelength. The model fitting with the closure amplitude 
constraints is free from any antenna-dependent amplitude errors, and decides on the best 
fitted model by searching over all the possible models. This proved highly effective in 

dealing with the Sgr A* observations when the radio emission of Sgr A* can be 
represented by an elliptical Gaussian of 3 parameters, i.e. sizes of major axis (0 major ) and 

minor axis (0 minor ) and, the position angle (PA) of the major axis (East of North). This is 

because that the closure phases, as a direct measure of the asymmetry in source structure, 



are measured to be consistent with a zero value. There are some data points deviating 
from the zero closure phases but having larger uncertainties. We thus restrict our model 
fitting to a single Gaussian component whose brightness distribution is symmetric with 
zero visibility phases. So, only the visibility amplitude is used in the model fitting. 
Before fitting, we performed a bias correction by solving two equations of the second and 

fourth moments of the measured visibility amplitude " . The best-fit model of the source 
structure is obtained by minimizing j 2 which is the weighted sum of the squared 

difference between the model and the measured (bias-corrected) visibility amplitudes 
over all the available baselines throughout the whole observation. In practice, a thorough 
search for minimum chi squares j 2 m is carried out in the 3-dimensional space that covers 
a wide range of the three model parameters (Q major , 0 minor and PA). This procedure is 

applicable to any source whose brightness distribution is characterized as a single 
Gaussian component. 

For the error estimation, in case that xl at X mm * s larger than unity, we scale up the 

68.3% confidence region (corresponding to the formal standard error lcr) as an increase 
of 7 2 from to j 2 in +Aj 2 withAj 2 = zL /N dof . Here, N dof is the summation of 

the difference between the number of visibilities N vis and the number of antennas N ant 

over all scans, minus the number of fitting parameters (3 for an elliptical Gaussian 
model), i.e. N dof = ^ ( N vis - N ant ) - 3 . By projecting this confidence contour onto the 

axis of parameter of interest, we can finally obtain lcr for that single parameter. 

Details of Revision of the Wavelength-dependent Scattering Law 

The consistent departure of the fitted 7 mm apparent major axis size from the pure 
scattering size predicted by the current scattering model seen in all seven epochs (Table 1) 
seems to imply that the intrinsic source size of Sgr A* at 7 mm has already added 
significantly in quadrature to the scattering size in producing the observed deviation. 



However, this conclusion explicitly assumes that the adopted wavelength-dependent 
scattering law is accurate enough. To check whether the difference is really due to 
intrinsic structure effect or an illusion caused by the inaccuracy of the current scattering 
model, we carry out weighted least-squares fit to the size measurements as a function of 
the observing wavelength. The weighting function is inversely proportional to the 
variance of the fitted size in Table 1. To minimize any temporal structural variability, all 
data points used except that at 3.5 mm are from near-simultaneous observations in 
February 1997 with the VLB A plus one antenna of the VLA. 

Supplementary Figure 1 and Supplementary Table 1 summarize various trials. The 
power-law form 0 = 0 lcm 7T was adopted in the fitting. The different trials made use of 

various subsets of the data in order to test the impact of including the shortest wavelength 
( A min ) as well as the longest wavelength (T max ) data, and also allowing the power law 

index a to be a free parameter in the fitting or fixed to the value of 2 as predicted by the 
scattering theory. For the major axis, we first tried to fit four subsets with A mm = 0.348, 
0.695, 1.350 and 1.953 cm, respectively, and T max =6.02 cm. The best-fit model of 
0major = (1.51 ± 0.06)T‘ 88±0 05 (errors are la) is obtained from the subset with A mm = 1.953 

cm . With a constant a = 2, two subsets with A mm = 1.350 and 1.953 cm, can be fit by 
the same 0 major = (1.38 + 0.02) A 2 . Thus, it is clear that the inclusion of 0.694 cm data did 

not lead to the best fitting in these fits for the major axis. The scattering size at 6.02 cm is 
heavily resolved by the VLBA. This severely limits the usable dataset. In particular, the 
fitted minor axis size could be arbitrarily small due to the poor constraint from the VLBA 
data along the North-South direction. As such, we tried the model fitting to a new sub¬ 
dataset with A- = 1.953 cm and /L., v =3.564 cm. The fitting results are 
® major = (1-39+ 0.05)A 200±a06 (varying a ) and 0 major = (1.39 ± 0.02)T 2 ( a = 2 ), 

respectively. Both fitting results are much better than above-mentioned fits. And, the 
fitted amplitudes 0 lcm from both fits are identical and, the fitted exponential index 



a = 2.00 + 0.06 is indistinguishable from the fixed number of 2. We thus conclude from 
these tests that 0 major = (1.39 ± 0.02)A 2 is the best-fit scattering law along the major axis 

direction. By assuming that the l 2 -dependence derived for the major axis also applies to 

the minor axis, we can obtain the best-fit wavelength-dependent minor axis size 
®minor = (0.69 ± 0.06)/l 2 from fitting to the subset of data with /L min = 1.35 cm. This 

newly revised wavelength-dependent two-dimensional scattering structure is also 
consistent with results from fitting to the closure amplitudes of somewhat different 
databases (STable 3 in ref. S4). 
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Supplementary Table 1. Summary of power-law fits to the major axial size of Sgr A* as 
a function of the observing wavelength. 
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Although it is widely accepted that most galaxies have supermassive black holes 
(SMBHs) at their centers 13 , concrete proof has proved elusive. Sagittarius A* (Sgr 
A*) 4 , an extremely compact radio source at the center of our Galaxy, is the best 
candidate for proof 6 ' 7 , because it is the closest. Previous Very Long Baseline 
Interferometry (VLBI) observations (at 7mm) have detected that Sgr A* is ~2 
astronomical unit (AU) in size 8 , but this is still larger than the “shadow” (a 
remarkably dim inner region encircled by a bright ring) arising from general 
relativistic effects near the event horizon 9 . Moreover, the measured size is 
wavelength dependent 10 . Here we report a radio image of Sgr A* at a wavelength 
of 3.5mm, demonstrating that its size is ~1 AU. When combined with the lower 
limit on its mass 11 , the lower limit on the mass density is 6.5xl0 21 M sun pc' 3 , which 
provides the most stringent evidence to date that Sgr A* is an SMBH. The power- 
law relationship between wavelength and intrinsic size (size oc wavelength 109 ), 
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explicitly rules out explanations other than those emission models with stratified 
structure, which predict a smaller emitting region observed at a shorter radio 
wavelength. 

Past VLBI observations 1216 of Sgr A* have revealed an east-west elongated 

structure whose apparent angular size at longer wavelengths is dominated by the 
interstellar scattering angle, i.e. 0 obs = ©^"'z 2 , where A is the wavelength in cm, and 

0 obs is the observed size in milli-arc second (mas). Thus, VLBI observations at shorter 

millimetre wavelengths where the intrinsic structure of Sgr A* could become 
comparable to the pure scattering size, are expected to show deviations of the observed 
size from the scattering law. This has been demonstrated by the recent detection of the 
intrinsic size at 7 mm (ref. 8). On November 3, 2002, we successfully carried out an 
observation of Sgr A* with the Very Long Baseline Array (VLBA) at its shortest 
wavelength of 3.5 mm (ref. 10). Our observation, with the steadily improved 
performance of the VLBA system, has produced the first high-resolution image of Sgr 
A* ever made at 3.5 mm (Fig. 1), which exhibits an elongated structure too. 

To yield a quantitative description of the observed structure, we tried a model 
fitting procedure 17 in which the amplitude closure relation is applied. Compared to the 
conventional VLBI self-calibration and imaging technique, this can improve on the 
calibration of current VLBI observations of Sgr A* (see Supplementary Information). 
We have applied this model fitting procedure to twelve sets of VLBA observations of 
Sgr A* made at a variety of wavelengths from 6 cm to 3.5 mm over the time range from 
1994 to 2004. Two experiments at 7 mm in 1994 are from VLBA archive. Table 1 lists 
the fitting results. A consistent position angle (-80°) of the scatter-broadened image can 
be seen in all the datasets, regardless of the observing epoch and wavelength. 
Furthermore, it is quite significant that in all seven experiments made at 7 mm, the fitted 
apparent major axis size is always larger than the largest known scattering size of 0.69 
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mas extrapolated from the existing scattering models. We thus perform weighted least- 
squares fit to the near-simultaneous angular size measurements (in February 1997) as a 
function of the observing wavelength (see Supplementary Information for more details). 

We conclude that the best-fit two-dimensional scattering structure is 
©maj 0r = (1 -39 ± 0.02)A 2 by 0 minor = (0.69±0.06)/l 2 with a position angle of -80°. This 

gives an even smaller scattering angle along the major axis direction. An immediate 
important conclusion is that the discrepancy seen in all 7-epoch 7 mm VLBI 
observations is real, implying the appearance of the intrinsic source structure at 
wavelengths of 7 mm and shorter (see Supplementary Figure 1). 

At 7 mm, two measurements from the March 2004 observations are used to get 
averaged sizes of the major and minor axes of 0.724 ±0.001 and 0.384 ±0.013 mas, 
respectively with position angle 80.6^ degrees. The difference between the measured 

and the extrapolated scattering sizes along the major axis is A0 = 0.053 ± 0.010 mas, 
significant at > 5cr level. By subtracting in quadrature the scattering angle, this suggests 
an intrinsic size of 0.268 ± 0.025 mas for the major axis. Similarly, we can derive an 
intrinsic size for the minor axis to be 0.190 ±0.057 mas, comparable to the major axis 
size. However, it should be kept in mind that the deviation seen for the minor axis is 
only significant at 1.6cr level. Note that the derived source size has greater statistical 
significance than the deviation, because the scattering size has been deduced with good 
accuracy. 

At 3.5 mm, the fitted apparent source structure from November 2002 VLBA 
observations is 0.21 'y® mas by 0.13 mas with a position angle of 79° (see 

Table 1). Thus, for the first time an intrinsic size of 0.126 ±0.017 mas for the major 
axis can be obtained from the deviation A0 = 0.042 ± 0.010 mas (at > 4cr level) in the 
measured major axis size from that of scattering angle at 3.5 mm. Mainly due to the 
limited resolution, the minor axis measurement at 3.5 mm, however, is inadequate to 
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make any firm claim on the determination of its intrinsic size. So, we defer any estimate 
of intrinsic minor axis size at 3.5 mm for future investigation. Past 3.5 mm VLBI 
observations with the heterogeneous Coordinated Millimeter VLBI Array, severely 
limited by its low sensitivity, could not warrant a model more complex than the circular 
one 18 . The best-fit circular Gaussian has a diameter of 0.18 ±0.02 mas, which is 
indistinguishable from the scattering size along the major axis and cannot give a 

Q 

meaningful estimate of the intrinsic structure . 

Thus we have sampled a zone of the SMBH closer to the event horizon than ever 
before, by detecting the intrinsic size of Sgr A* to be only 1.01 AU at a distance of 8.0 
kpc, or 12.6R SC , where R sc (= 1.2 x 10 12 cm) is the Schwarzschild radius of a 4x 10 6 M sun 

SMBH. By assuming a spherical structure, we obtain a lower limit to the mass density 
of Sgr A* of 6.5 x 10" M sun pc' . Here, we used the lower bound to the mass of Sgr A*, 
derived from the upper limit to the intrinsic proper motion of Sgr A* itself 11 , which is 
about 10% of the 4xl0 6 M sun inferred from the stellar orbital motions 6 ’ 7 , and the upper 
limit to the source intrinsic size (this work). This mass density is at least 4 orders of 
magnitude greater than that determined from dynamical measurements of stellar 
velocities 6 . This is because here we are probing directly the structure of Sgr A*, where 
the assumed mass estimate refers to the value within Sgr A*. We note also that this 
mass density is almost 12 orders of magnitude greater than the estimate for NGC 4258 
(ref. 19), one of the best known SMBHs. Such an extraordinarily high mass density 
robustly rules out the possibility of Sgr A* being a compact dark cluster of stellar 
remnants as it would have an unreasonably short lifetime of less than 100 yrs (ref. 20), 
and thus argues strongly in favor of the SMBH nature of Sgr A*. To prove that Sgr A* 
is indeed an SMBH requires an unambiguous demonstration that Sgr A* possesses an 
event horizon. It is intriguing that the detected intrinsic size at 3.5 mm is about two 
times the diameter of the shadow caused by the strong gravitational bending of light 
rays 9 . Thus, it is very promising that VLBI observations of Sgr A* at 1 mm or shorter 
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will reach the region comparable to its shadow, which can be used to differentiate 

9 21 

between the SMBH scenario and other supermassive non-baryonic stars ’ . 

The two-point fit to the well determined intrinsic sizes at 7 and 3.5 mm shows a 
A p -dependence of the intrinsic source size with ft = 1.09 ^ (Fig. 2). Also plotted in 

Fig. 2 are two inferred lower limits of 0.02 and 0.008 mas to the intrinsic size at 1.3 and 
0.8 mm, respectively, from the absence of refractive scintillation' . These lower limits 
are consistent with the extrapolation of the A? -dependence. However, we note that 
these two lower limits are only about 2 and 0.8 R sc , which are smaller than the last 
stable orbit (LSO) radius of 3 R sc for a non-rotating (Schwarzschild) black hole. There 

IT 

is some evidence that Sgr A* is a rotating black hole' . For a prograde maximally 
rotating Kerr black hole the LSO radius is 0.5 R sc . The LSO establishes the lower limit 
to the emission region size. Hence the A? -dependence will eventually reach a minimum. 
As such, the turn-over frequency' seen in the entire spectrum of Sgr A* might tell us 
the smallest size of the emission, which can be further used to constrain its spin, if any. 

The extrapolated intrinsic size at 1.35 cm is 0.555^'|^ mas. This, when compared 
to the scattering angle of 2.576 ±0.036 mas, is consistent with the idea that the 

scattering effect dominates the observed source size at 1.35 cm. This deduced source 
size is also formally consistent with the reported detection of 0.726 ]',j mas, to within 

o 

the uncertainties . 


The derived A 109 -dependence requires that the emission at different wavelengths 
is dominated by different emitting region and thus conclusively exclude those models 
without the stratified emission structure. Along with the detected intrinsic major axis 


size, we can derive a lower limit to the intrinsic brightness temperature as 
, K; here S A is the flux density in Jy at wavelength A in cm and 


T b > 1.36x10" x 


su 2 


0 


0 int is the intrinsic major axis size in mas. There is a wavelength dependence of the 
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lower limit T b as A ° 018 (assumingS^ oc A “ ). Using the flux densities of 1.0 and 1.2 Jy 
at 7 and 3.5 mm, respectively, the corresponding minimal T b is 0.9xlO 10 and 1.2 xlO 10 

K, greater than the prediction of the spherical accretion model . However, this lower 
limit of 10 10 K, and the spatial distribution of the radio emission, can be explained 
easily by the inhomogeneous jet model 26 ' 27 , in which the magnetic field and the electron 
number density vary with the distance (r) to the origin of the jet as r' 1 and r' 2 , 
respectively. On the other hand, the radiatively inefficient accretion flow (RIAF) model 
of Sgr A* (ref. 28) can also account for a brightness temperature of >10 10 K as well as 
the observed spectral energy distribution. The prediction of A 0 ' 9 from the hybrid 
thermal-nonthermal synchrotron radiation from RIAF (ref. 29) is in agreement with the 
estimated A 032 relation. Here, the possible existence of strong outflows from the 
accretion disk was not taken into account. To further discern between them, it is 
important to study the correlation between the detected X-ray variability and the 
variations frequently seen in the radio to sub-millimetre wavelengths, which would 
yield further information on the intrinsic density structure of the emitting zone. 
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Figure 1: The first high-resolution VLBI image of Sgr A* at 3.5 mm obtained with 
the VLBA on November 3, 2002. The observations were dynamic scheduled to ensure 
the good weather conditions at most sites, and the data were recorded at the highest 
possible recording rate of 512 Mbps (Mega Bits Per Second). Standard visibility 
amplitude calibration including the elevation-dependent opacity correction was done, 
and the final image was obtained after several iterations of the self-calibration and 
cleaning procedures. The calibrated total flux density is about 1.2 Jy. a, A uniformly 
weighted image with the restoring beam (indicated at the lower left comer) of 1.13 mas 
x 0.32 mas at 9°. The peak flux density is 1.08 Jy beam" 1 . Contour levels are drawn at 
3cr x (-1, 1, 2, 4, 8, 16, 32), 3cr=17.5 mJy beam" 1 , b, A super-resolution image with a 
circular beam of 0.20 mas from which an east-west elongated structure can be seen (see 
Table 1). Note the different scales. The contour levels are the same as that in a with the 
corresponding peak flux density of 1.01 Jy beam' 1 . 

Figure 2: Intrinsic major axis size vs. observing wavelength. The solid line represents 
the two-point fit from the detected intrinsic sizes at both 3.5 and 7 mm (this work). Also 
plotted are the lower limits to the intrinsic sizes at 1.3 and 0.8 mm (ref. 22), and the 
reported detection at 1.35 cm (ref. 8). The extrapolated intrinsic sizes at 1.3 and 0.8 mm 
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are about 4.3R SC and 2.5R SC , respectively. Here, R sc (= 1.2 x 10 12 cm) is the 
Schwarzschild radius for a 4x 10 6 M sun SMBH. 

Table 1. Parameters of the elliptical Gaussian model for Sgr A* 


2 

v (bw, rate ) 3 

Epoch 

0 . M 
major 

0 . b ’ d 
minor 

PA d 

(cm) 

(GHz) (MHz, bits) 

(dmy) 

(mas) 

(mas) 

(°) 

0.348 

86.236(128,2) 

03Nov2002 

0 21 +0 02 
U - Z ' A -0.01 

0 13 +005 

u * 1 J - 0.13 

79 !; 2 

0.694 

43.213 (32,2) c 

14Febl997 

0 7l +0 ' 01 

n 49 +0 - 05 

74 ! 2 

0.695 

43.175 (128,2) 

08Mar2004 

O 999 + 0-002 

U. 1 QQ02 

0 395 +0 019 

80.4™ 

0.695 

43.175 (128,2) 

20Mar2004 

rv '79^+0.002 
U. /Z 9 _o.o02 

O 999 +0.020 
U .9 / _o.oi8 

80.8™ 

0.695 

43.151 (64,1) 

26Aprl994 

0 72 +001 

0 39 +0 07 

y - 0.07 

78! 2 

0.695 

43.151 (64,1) 

29Sepl994 

0 72 +001 
' ^-0.01 

0 42 +0 03 

^•^- 0.03 

79 

0.695 

43.135 (32,2) 

24Aprl999 

0 69 +0 01 

0 33 +0 04 

u,JJ - 0.04 

83 

0.695 

43.135 (32,2) 

23Mayl999 

0 71 +0 01 
' A -0.01 

0 44 +0 - 02 
u *^ + -0.02 

19 + _\ 

1.350 

22.229 (32,1) c 

12Febl997 

r \ ro + 0.06 

Z.J^_o .05 

1 4 <+ 0-23 

A ‘^ - 0.38 

S 3 U 5 

1.953 

15.361 (32,1) c 

12Febl997 

J. 9 9 _oo7 

2 . 70 ™° 

83! 3 3 

3.564 

8.417 (32,1) c 

07Febl997 

17 5 +0,5 

8 . 51 ^ 

87 ! 3 

6.020 

4.983 (32,1) c 

07Febl997 

43.01“ 




Results from the model fitting procedure which implicitly uses the amplitude closure relation (see 
Supplementary Information for details). Also listed are some details of the VLBI observations. Except 
five observations in February 1997 which used the VLB A and one VLA antenna, all the other 
observations were performed by the VLBA. 
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a the observing frequency in GHz. Numbers in parenthesis are the recording bandwidth in MHz and the 
sampling rate (1 bit or 2 bits). With the same total bandwidth, the noise level from 2-bit quantization 
data is about V2 times lower than that from 1-bit quantization. The recording rate is two times the 
product of the bandwidth and the sampling rate. 

b full width half maximum (FMHM) 

c dual polarization observation. Both left and right circular polarization data have the same recording 
mode (listed) and, both data were used in the model fitting. 

d in all cases, errors are 1 a 



